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Abstract. We show strong evidence for the absence of a finite-temperature spin 
glass transition for the random-bond Ising model on self-dual lattices. The analysis is 
performed by an application of duality relations, which enables us to derive a precise 
but approximate location of the multicritical point on the Nishimori line. This method 
can be systematically improved to presumably give the exact result asymptotically. 
The duality analysis, in conjunction with the relationship between the multicritical 
point and the spin glass transition point for the symmetric distribution function of 
randomness, leads to the conclusion of the absence of a finite-temperature spin glass 
transition for the case of symmetric distribution. The result is applicable to the random 
bond Ising model with ±J or Gaussian distribution and the Potts gauge glass on 
the square, triangular and hexagonal lattices as well as the random three-body Ising 
model on the triangular and the Union- Jack lattices and the four dimensional random 
plaquette gauge model. This conclusion is exact provided that the replica method is 
valid and the asymptotic limit of the duality analysis yields the exact location of the 
multicritical point. 
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1. Introduction 

Properties of finite-dimensional spin glasses are still under active current investigations 
after thirty years since mean field analyses of the basic model [Il|2]. The difficult problem 
of whether or not the mean-field predictions apply to realistic finite-dimensional systems 
is still largely unsolved. One of the outstanding problems is the existence or absence of 
the spin glass phase. Most of the current investigations on this problem are carried out 
using numerical methods [3, SI [5]. 

Two typical spin glass models have been examined extensively, the Gaussian and 
±J Ising models. In three dimensions, researchers have arrived at the consensus that 
there are finite-temperature spin glass transitions in both models. On the other hand, 
in two dimensions, numerical investigations show evidence that there would be no 
finite-temperature spin glass transition in both models P, [3, El El [101 [HI [12], [13], [H]. 
Unfortunately no reliable analytical evidence for the problem of the existence or absence 
of the spin glass phase in finite dimensions has been established. 

Very little systematic analytical work for finite-dimensional spin glasses exists. An 
exception is a technique based on the gauge symmetry to derive the exact value of the 
internal energy, a rigorous upper bound of the specific heat, several set of rigorous 
inequalities and exact relations in a special subspace, known as the Nishimori line 
[T5l [16] . In the present study, we develop an argument by the gauge symmetry in 
conjunction with the duality and the replica method to study the problem whether or 
not a finite-temperature spin glass transition exists in two dimensions. We analyze the 
problem of the spin glass transition point for the ±J Ising model and the Gaussian 
Ising model on self-dual lattices by means of the duality [T71 [18]. The theory is 
applicable directly to the square lattice, but the triangular and hexagonal lattices can 
also be reduced to be self-dual using the duality in conjunction with the star-triangle 
transformation [T9|. In the present study, we arrive at the conclusion that no finite- 
temperature spin glass transition exits in the symmetric distribution of randomness, for 
example, p = 1/2 for the ±J Ising model. The result is justified under the validity of 
the replica method. It should also be remembered that the prediction of the duality 
method for the transition point is expected to be exact only in the asymptotic limit of 
large cluster size in the sense we shall define in the following sections. 

This paper is organized as follows. In the next section, we recall the basic 
formulations of the duality in spin glasses with the replica method and the relation 
derived by the gauge symmetry to set a stage to show the absence of a spin glass phase 
in the following section. We develop an analytical argument to derive our result in 
section 3. In this section, we show the limit of applicability of the present result. The 
final section is devoted to conclusion and discussions. 
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Figure 1. Self-dual lattices. (A) Square lattice. (B) Triangular lattice. (C) 
Hexagonal lattice. (D) Hierarchcial lattice. The triangular and hexagonal 
lattices become self-dual by means of the duality combined with the star- 
triangle transformation |19j . Construction of a hierarchical lattice starts from 
a single bond, and we iterate the process to substitute the single bond with the 
unit cell of a complex structure [20l [211 112] • 



2. Duality and gauge symmetry 

Let us review several known facts in the present section to fix the notation and prepare 
for tlie developments in the next section. The following arguments are applicable to 
several spin glass models on self-dual lattices as shown in figure [H We take the ± J 
Ising model on the square lattice as an example here for simplicity. It is straightforward 
to generalize the following arguments to other spin glass models as will be explained in 
the next section. 



2.1. ±J Ising model 

The Hamiltonian of the ±J Ising model is given by 

H = -"^JijSiSj, (1) 

(ij) 

where Si is the Ising spin taking ±1, and Jij{= ±J) denotes the quenched random 
coupling. The summation is taken over nearest neighbouring pairs of sites. The partition 
function with fixed randomness is 

ZiK,{T,,}) = J2I[e^^'^'''^, (2) 

{s^} (ij) 

where K = (3 J is the coupling constant, and Tjj(= ±1) is the sign of the random coupling 
Jij. The distribution function of Tij is 

Pir.^) = V5{r,j - 1) + (1 - V)5{n, + 1) = ^fj^K, ' 
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where we defined Kp by exp {—2Kp) = (1 — p)/p. 
2.2. Duality and multicritical point of spin glasses 

We here review the analysis of the location of the multicritical points by the duality 
in spin glasses [23l [2^ . The multicritical point is connected with a critical point for 
p = 1/2 {Kp = 0) by a relation introduced in section 2.4. 

We apply the replica method to the ±J Ising model. The n-replicated partition 
function after the configurational average is 

Zn{K. Kp) = 2k ^ n e^^^-Z{K, {r.,})^ (4) 

where n stands for the number of replicas. We generalize the duality argument to 
this n-replicated ±J Ising model [23l [2lj. For this purpose it is useful to define the 
edge Boltzmann factor Xk {k = 0,1, ■ ■ ■ ,n), which represents the configuration-averaged 
Boltzmann factor for interacting spins with k antiparallel spin pairs among n nearest- 
neighbour pairs for a bond (edge). The duality gives the relationship of the partition 
functions with different values of the edge Boltzmann factor as given by 

The dual edge Boltzmann factors xl are defined by the discrete multiple Fourier 
transforms of the original edge Boltzmann factors, which are simple combinations of plus 
and minus of the original Boltzmann factors in the case of Ising spins. Two principal 
Boltzmann factors Xq and Xq are the most important elements of the theory 

cosh {Kp + riK) 
cosh K, 



MK.K)^ -"::r-\ (6) 



xl{K,Kp) = [V2 cosh kY . (7) 

We extract these principal Boltzmann factors from the partition functions in equation 
(IS]) to measure the energy from the all-parallel spin configuration. The result is written 
as, using the normalized edge Boltzmann factors Uj = Xj/xo and u* = x*/xq, 

Xq{K, Kp)^''Zn{Ui,U2, ■■■ ,Un) = Xq{K, Kp)^'' Zn{u\,ul, " " " , <), (8) 

where Zn{ui, ■ ■ ■) and Zn{u{, ■ ■ ■) are defined as Zn/x^^ and Zn/ {xq)^^ . 

The duality identifies the critical point under the assumption of a unique phase 
transition. The critical point is given as the fixed point of the duality transformation 
and is known to yield the exact critical point for a simple ferromagnetic system on 
the square lattice [T71 [18]. In order to obtain the multicritical point of the present 
replicated spin glass system, we set K = Kp, which defines the Nishimori line (NL) on 
which the multicritical point is expected to lie. Since IS Si multivariable function, 
there is no fixed point of the duality relation in the strict sense which satisfies n 
conditions simultaneously, ui{K) = ul{K),U2{K) = U2{K) , ■ ■ ■ , Un{K) = u^{K). This 
is in sharp contrast to the non-random Ising model, in which the duality is a relation 
between single-variable functions. We nevertheless set a hypothesis that a single 
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cluster 1 cluster 2 



Figure 2. A cluster is the unit plaquette encircled by white spins. The spins marked 
black on the original lattice are traced out to yield interactions among white spins in 
clusters. 



equation xo{K,K) = Xq{K,K) gives the location of the multicritical point for any 
replica number n [231 l2ll [251 ISSl [271 12H] • The quenched limit n — > for the equation 
xo{K, K) = x*q{K, K) then yields 

-plogp- (1 -p)log(l -p) = ^log2. (9) 

The solution to this equation is Pc = 0.889972(^ 0.8900). 

This value of Pc is very close to numerical results but shows small deviations on 
several self-dual hierarchical lattices, for which numerically exact results can be derived, 
as in tabled] [2H]- Such discrepancies may come from the following two facts. One is 
that the condition Xq{K, K) = Xq{K, K) is different from the strict fixed-point condition 
and the other is that we have considered only the quantity defined on a single bond, 
the principal Boltzmann factors Xq and Xg, which does not necessarily reflect the effects 
of frustration inherent in spin glasses. The improvements with these points taken into 
account will be explained in the next section, following [291 EO] . 

2.3. Improved method 

As shown in figure El let us consider to sum over a part of the spins, called a cluster 
below, on the square lattice to deal with the effects of frustration rather than a single 
bond considered in the naive approach described in the previous section. Then, the set 
of clusters must be chosen to cover the whole lattice under consideration as in figure [21 
where two examples of the choice of clusters are depicted. The partition function can 
now be expressed as a function of the configuration of spins around a cluster (marked 
white in figure [2]) . The duality is written as 

^i^H4^\4^---) = 4^n4^^\a;^(^\---)- (10) 

The superscript s stands for the type of the cluster that one chooses. The quantity 
x^'*^ is the local Boltzmann factor including many-body interactions generated by the 
summation over spins marked black in figure [21 We define the principal Boltzmann 
factors Xq^ and its dual those with all spins surrounding the cluster in the up 
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state. We assume that a single equation gives the accurate location of the multicritical 
point for any number of n, similarly to the naive conjecture, 

xi'\K,K)=xt\K,K). (11) 

This is the improved method to predict the location of the multicritical point with higher 
precision than the naive conjecture. 



Pc (conjecture) Pc (improved) Pc (numerical) 



0.8900 
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B920 


O.i 


B915(6) 


0.8900 


O.J 


B903 


O.i 


B903(2) 


0.8900 


O.i 


B892 


O.i 


B892(6) 


0.8900 


O.J 


B895 


O.i 


B895(6) 


0.8900 


0.1 


B891 


O.i 


B890(6) 



Table 1. Comparison of the naive conjecture, improved method and numerical 
estimations for self-dual hierarchical lattices [29] , 

The improved method for the ±J Ising model indeed has given the results in 
excellent agreement with the exact estimations within numerical error bars on several 
self-dual hierarchical lattices as summarized in table [1] [29l [30]. In addition, recent 
numerical investigations on the square lattice have given pc = 0.89081(7) [31] and 
Pc = 0.89061(6) [32], while the improved method has estimated Pc = 0.890725 by cluster 
1 of figure E], and pc = 0.890822 by cluster 2 [30]. If we deal with clusters of larger sizes, 
the improved method should show systematic improvements toward the exact answer 
on the location of the multicritical point. We use this property of the improved method 
to show that the spin-glass transition temperature would be zero on self-dual lattices, 
which means the absence of a finite-temperature spin glass transition. 

2.4- Relation between different replica numbers 

Another important piece of information is the relation that the (n-|- l)-replicated system 
with p = 1/2 {Kp = 0) is equivalent to the n-replicated system on the NL [33] . 

After the gauge transformation and summation over gauge variables, the 
exponential in equation (jl]) turns to a partition function with coupling Kp [151 [IS]: 

K.) = 2N.(2JshK)^s ^ ^"--^^ ■ ^^^^ 

where Ng is the total number of spins. It readily follows from this equation that 
Zn{K,K) and Zn+i{K,0) are essentially equal to each other: 

2^' {2 cosh K)^^Zn{K,K) = 2^^+^^Z„+i(ir, 0). (13) 

This result is quite general since we did not use the properties of a specific lattice. 
If there is a singularity in the partition function Zn{K,K) on the left-hand side at 
some point [Kc, Kc) for replica number n, Zn+i{K, 0) on the right-hand side has also a 
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singularity at {Kc, 0) for n + 1. According to this relation, it is sufficient to evaluate the 
location of the multicritical point in the limit — 1 in order to study whether or not 
there is a finite-temperature spin glass transition for Kp = in the quenched disorder 
system (n — > 0). In the following section, we will derive clear evidence for the absence of 
a finite-temperature spin glass transition for the case of p = 1/2 {Kp = 0) by estimating 
the location of the multicritical point in the limit n — 1 by the improved method. 
The present argument is valid as long as the replica method is reliable, in particular in 
the limit n — — 1. 



3. Absence of a finite-temperature transition for p = 1/2 

In this section we show that Tsg = for p = 1/2 {Kp = 0) for an arbitrary cluster of the 
improved method introduced in the previous section. Since the improved method would 
give the exact result in the limit of infinitely large clusters because we take the full trace 
over all spins in the system, except for the boundary spins, to give the exact partition 
function, we expect our conclusion Tsg = to be valid not as an approximation but as 
the exact conclusion. 

According to the improved method, the location of the multicritical point for the 
n-replicated ±J Ising model is given by equation (ITT]) . The explicit expressions of two 
principal Boltzmann factors are given as 

4'\k, Kp) = \ E n e^^^^^ ( E ' n '^^'^'''^ I , (14) 

{2coshKp)^B ^^^^^.^ [^^^^ ^^^^ J 



and 



{2cos\iKpYB ^^^^^.^ 



y^^{S^} (ij) 




(15) 

-(^^^^£n'e-^ (16) 

where the prime on the summation denotes the condition that all the spins at the 
perimeter of the cluster are up and the prime on the product represents that the pairs 
are restricted to those within the cluster. The quantity A'^^^' is the number of bonds in 
the cluster. 

The principal Boltzmann factor Xq^ is regarded as the n-replicated partition 
function of the finite-size cluster after the configurational average, which is cut out 
of the self-dual lattice under consideration: 

4'\K,Kp) = \ E I['^''^^''Z^'\K,{r,,}r. (17) 

(2coshKp)^B ^^^^^^^^ 
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Figure 3. Dual pairs for cluster 1 and cluster 2. The same symbols are used 
as in figure [2j 



Gauge transformation and summation over gauge variables reduce the exponential factor 
in this expression to a partition function of the cluster with coupling Kp [T5l [T6] , 



4-'(A',A'p 



where A^^'' denotes the number of spins inside the cluster under consideration. 



The quantity Z**^*-' in the dual principal Boltzmann factor Xq , in contrast, is not 
gauge invariant as seen in equation (fT5l) . We here apply the duality transformation to 
the cluster with all edge spins being up. We can then obtain another expression 
of as EU] 



dual 



{S^} (tj) 

■,Nf--N'=J-/2-l yis) 



(19) 



2-r-^^Br^-^Z^>{K,{n,}). (20) 

The duality in the cluster allows us to rewrite Z*^^\K, {Tij}) as the partition function 
Zj^\k, {Tij}) defined on the dual lattice of the cluster under consideration, denoted by 
the subscript D, as in figure [3l It is noted that the total number of spins in the dual 
cluster is given by the number of plaquettes in the original cluster, N^-. Again the 
technique of gauge transformation with N^- gauge variables is applicable to rewrite the 
dual principal Boltzmann factor 



is) 



iK,Kp) 



2n(N\ 



cl._^cl./2-l) 



-E4^)(ir„{r,,})4^^(K,K})". (21) 



2^^'(2coshirJ^-^' 



Our task is thus to solve the following equation derived from x, 
condition K 



is) 




2n{Nf--Nf/2-l) 



{Tij} 



Ez(^)(ir,{r,,}) 

{Tij} 



X, 



n+1 



*(s) 



with the 



(22) 



This equation gives the location of the multicritical point, which corresponds to the 
singularity on the left-hand side of equation ( fT3l) . To show the absence of a finite- 
temperature transition for Kp = in the quenched system (n 0), we consider the 
limit of n — > — 1 in the above equation, since the multicritical point for n — > — 1 is 
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Figure 4. Behaviours of Tsg as a function of the replica number n derived by 
Xg''^ = Xq*"^^ . The dashed curve denotes the results by the naive conjecture, and 
the solid curve represents those by the improved method for cluster 1. The 
curve showing the points given by cluster 2 coincides with that by cluster 1 in 
this scale. 

equivalent to the critical point for = in n — > due to the relation ( |T3ll . It should 
be noted that the two partition functions Z^^^ and do not have any singularity since 
they are given by the summation over spins of finite-size systems. 

Let us assume that there is a finite-temperature transition Tgc = l/-ft"sG- Then 
the partition functions Z^^'^ and z'^ on the original and dual clusters have some finite 
values. Equation fl22|) then reduces to, in the limit of n — > —1, 

2<-^sV2-i ^ ^ (23) 

The fact that A^'^'- ^ N^- /2 + 1 for any finite-size cluster cut out of the self-dual lattices 
leads us to the conclusion that equation (123|1 can not be satisfied. In this sense, the spin 
glass transition point of the ± J Ising model on a self-dual lattice goes to zero Tsg 
as n ^ — 1. By taking the asymptotic limit of large clusters, we expect the result to be 
exact. 

We show the behaviour of Tgc as a function of n derived from Xq — Xq 111 figure 
m for the naive conjecture and the improved methods. It is clearly observed that the 
value of Tsg depends on the degree of improvement for > — 1 but it is fixed to Tsq at 
n = —1. 

The above arguments are applicable as long as equation fl22|) can be established. In 
this formulation, we used the replica method and the self-duality of lattices. The present 
results are acceptable under the validity of the replica method and are applicable to 
self-dual lattices, the square and several self-dual hierarchical lattices. Using the star- 
triangle transformation, we can derive essentially the same relation as equation (l22l) for 
the triangular and hexagonal lattices. 

We can apply the above formulation to the Gaussian Ising model for Jq = 0, where 
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Jo denotes the mean of the distribution of randomness because essentially the same 
improved method works [30] and the same relation as equation ( IT3|) holds [33]. The 
same is true for the g-state Potts gauge glass and the random Ising model with three- 
body interactions on the triangular and the Union- Jack lattices. The multicritical point 
for these models can be also given by the improved method [231 EH EOj [36l [37] . For 
instance, the coefficient on the left-hand side of equation f|T3|) is slightly modified and 
written by the state number q instead of 2 for the case of the Potts gauge glass. Also, 
the duality structure of the four- dimensional random plaquette gauge model is exactly 
the same as the ± J Ising model on the square lattice, and therefore the present analysis 
applies without change [25l [26] . 

4. Conclusion and discussions 

We showed the absence of a finite-temperature spin glass transition for several spin 
glass models with symmetric distribution of randomness in finite dimensions by using 
the duality and gauge symmetry. This result is applicable to spin glass models with the 
Nishimori line, the ±J Ising model, the Gaussian Ising model, and the g-state Potts 
gauge glass on several self-dual lattices. We remark that the triangular and hexagonal 
lattices are included as the self-dual lattices by means of the duality in conjunction with 
the star-triangle transformation. The random Ising model with three-body interactions 
on the triangular and the Union- Jack lattices and and the four-dimensional random 
plaquette gauge model also have no finite-temperature spin glass transition for the 
symmetric distribution of randomness. 

Many researchers believe the absence of a finite-temperature spin glass transition 
for the random bond Ising model in two dimensions from numerical investigations. The 
present analysis lays an analytical foundation of this expectation. Although a real-space 
renormalization group calculation suggests a finite spin glass transition temperature for 
the triangular lattice [38] , we believe that the approximation involved there is too crude 
to be qualitatively reliable. 

We showed Tsg = for the symmetric distribution of randomness, considering 
an arbitrary degree of the improved method. The conclusion derived in the present 
analysis is exact under the validity of the improved method in the asymptotic limit of 
large clusters and the replica method. 

Our result is significant in the sense that this is the first analytical and systematic 
evidence for the conclusion of the absence of a finite-temperature spin glass transition 
of the random-bond Ising model in two dimensions and related systems with symmetric 
distribution function of randomness. We believe that further generalizations to other 
systems are worth the efforts. 
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